function OpenCRTBP_MP_BaseFigure
%% Constants
G           = 6.674e-11 * ((1/1000)^3);
mass_Mars   = 6.39e23;
mass_Deimos = 1.8e15;
mass_Phobos = 10.8e15;
GM_Mars     = G*mass_Mars;
GM_Deimos   = G*mass_Deimos;
GM_Phobos   = G*mass_Phobos;
GM_sun      = 1.327124400*10^(11);     
%[km^3/s^2] Gravitational Parameters

u_sm = GM_Mars / (GM_Mars+GM_sun);
u_md = GM_Deimos / (GM_Deimos+GM_Mars);
u_mp = GM_Phobos / (GM_Phobos+GM_Mars);
%[] Mu constant
%    sm --> Sun-Mars system
%    em --> Earth-Moon system

% DU_sm = 235.34e6; 
DU_mp = (9236.53+9517.58)/2;
% DU_md = (23455.5+23470.9)/2;
%[km] Distance from Sun to Mars
%[s] Time unit 
%[] Mu constant
%    se --> Sun-Earth system
%    em --> Earth-Moon system
rp = 10/DU_mp;
rm = 3389/DU_mp;
figure
hold on
[L1,L2,L3,L4,L5] = librationPoints(u_mp);
Lse = [L1,L2,L3,L4,L5];
for ii = 1:5
    plot3(Lse(1,ii),Lse(2,ii),Lse(3,ii),'bo','MarkerFaceColor','b','MarkerSize',2)
end
[xs,ys,zs] = ellipsoid(-u_mp,0,0,rm,rm,rm,50);
fv = stlread('Phobos.stl');
fv.vertices = (fv.vertices)/10/DU_mp;
fv.vertices(:,1) = fv.vertices(:,1)+(1-u_mp);
patch(fv,'FaceAlpha',1,'EdgeColor','none');
surface(xs,ys,zs,'FaceColor',[255 165 0]/255,'EdgeColor','None');
xlim([-1.3,1.3]);
ylim([-1.3,1.3]);
axis equal
set(gca,'FontSize',20) 
end